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Abstract. Stanley's non-negativity theorem is at the heart of many of the 
results in Ehrhart theory. In this paper, we analyze the root behavior of general 
polynomials satisfying the conditions of Stanley's theorem and compare this 
to the known root behavior of Ehrhart polynomials. We provide a possible 
counterexample to a conjecture of the second author, M. Beck, J. De Loera, 
J. Pfeifle, and R. Stanley, and contribute some experimental data as well. 



Let P be a convex polytope in W 1 with vertices in Z" and affine span of dimension 
d. We will refer to such polytopes as lattice polytopes and to elements of Z" as lattice 
points. By a remarkable theorem due to E. Ehrhart, jS], the number of lattice points 
in the t th dilate of P, for non-negative integers t, is given by a polynomial in t of 
degree d called the Ehrhart polynomial of P. In this paper we will investigate some 
differences between the root behavior of Ehrhart polynomials for elements of an 
arbitrary collection of polytopes of dimension less than or equal to d and the root 
behavior of an arbitrary collection of polynomials of degree less than or equal to d 
satisfying a certain non-negativity condition. 



We will begin by reviewing basic facts about Ehrhart polynomials. It is well 
known, e.g. chapter 4 of JBJ, that for a polynomial / of degree d over the complex 



Given a lattice polytope P. we denote its Ehrhart polynomial by Lp(t) and let 



denote the Ehrhart series for P. There is a well-known relationship between this 
encoding of Lp(t) and the polynomial itself, namely that Lp(t) can be expressed 
as 



This is easily seen by expanding the rational function as a formal power series. 
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1. Ehrhart Theory 



numbers there exist complex values h* so that 
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Thus, Ehrp(x) encodes the change of coefficients for Lp(t) corresponding to the 
change from the standard monomial basis to the basis 

It turns out that representing g(t) in this way can be very profitable. One of the 
most important known results about Ehrhart polynomials is the following theorem 
due to R. Stanley, known as Stanley's non-negativity theorem. 

Theorem 1.1. (see [7] and If P is a d- dimensional lattice polytope, then 
(h* ,...,h* d )€(Z> ) d +\ 

The non-negativity theorem is at the heart of much of what is presently known 
about roots and coefficients of Ehrhart polynomials. However, not every polyno- 
mial with non-negative integer hj's is an Ehrhart polynomial, hence we make the 
following definition. 

Definition 1.2. A non-zero polynomial satisfying the condition that 

(h* ,...,h* d )e (K>o) d+1 

is called a Stanley non-negative, or SNN, polynomial. 

In this paper we will see some of the differences between the root behavior of 
Ehrhart polynomials and of arbitrary SNN polynomials. 

2. Norm Bounds and Growth Rates 

In this section we review a norm bound on roots of SNN polynomials and some 
results and conjectures about growth rates of roots of SNN and Ehrhart polynomi- 
als. In £Q, it was shown that for polytopes of fixed dimension d the roots of Lp(t) 
are bounded above in norm by 1 + (d+ 1)!. It was further suggested that this might 
be made polynomial in d. In response, the first author proved the following. 



Theorem 2.1. (see /// is an SNN polynomial, then all the roots of f lie inside 
the closed disc with center ^ and radius d{d — -|) . 

The proof of this can be found in 0], but we will find it useful to sketch the 
argument. If we represent such an / as 

/(*) = E ft 5 



t + d-j 
d 



3=0 

then to evaluate / at t S C we take a non-negative linear combination of the d+1 
points ( t+ j _ "') S C, < j < d. If all of these points arc in a common half space H 
of C with zero on the boundary, then f(t) ^ 0. Thus, one only needs to show that 
this is satisfied for t ^ {z : \z + i| < d(d - |)} to prove Theorem 12. II 

The above bound is essentially optimal, as the following theorem demonstrates. 

Theorem 2.2. (see j The polynomial Sd(t) = J2j=o C + d _J ) * s an Ehrhart poly- 
nomial whose roots all have real part = ^-. Further, if a d is the root of Sd(t) of 
maximal norm, then 

1 d{d + 2) 
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Thus, a norm bound for Ehrhart polynomial roots cannot be better than qua- 
dratic in d. It was suggested in [2] that Sd{t) possesses the roots of maximal norm 
among all dimension d polytopcs with interior lattice points, and this was proved 
for d = 2, 3. In response to the analogous question for SNN polynomials, wc offer 
the following theorem and conjecture. 

Theorem 2.3. For the polynomial Md(t) = (* "t ) + which is not an Ehrhart 
polynomial, if (3d is the root of M<j(i) of maximal norm, then 

\Pd+h = -+o(i), 



Proof. Note that the following closely follows the proof of Theorem 12 . 21 given in 

For any lattice polytope P, h* d is equal to the number of interior lattice points in 
P. If this is non-zero, then h*, which records Lp(t) — d+ 1 when P is d-dimensional, 
must also be non-zero. As this condition is not satisfied by Md{t), Md(t) is not an 
Ehrhart polynomial. 

By a result of Rodriguez- Villcgas in [5] , since the roots of the numerator of the 
generating function for Md(t) lie on the unit circle, all the roots of M d {t) are on 
the line x = = ^-. If s = ^ + hi, b > 0, is a root of M d (t), then we have 

(2.1) (s + d)(s + d - 1) • • • (s + 1) = -s(s - 1) • • • (s - d + 1), 

as any root s of Md {t) satisfies 

' s + d\ ( s 



d J \d 

Writing s — j = Sj = rje l8j and noting that |s +j + 1| = \s — j\ implies s + j + 1 = 
rje l ^~ 6 i\ we can rewrite H2.1|l as 

r_l\d+i _ e i(2e +---+2e d _i) _ 

We now substitute \ + 4>j = 6j, where <fij £ (0, ^]. This gives a new equation, 
Therefore, we must have, for some positive odd value of k, 

, d-l 
Z 

By definition, cot 0, = -rj-j- for j = 0, ...,d— 1. Thus, s is a root of Md{i) of 
maximal imaginary part if and only if 

d-i 

(2.2) p d (b) := ^ cot- 

3=0 

as each pd(b) is a strictly decreasing function of b. Say that Pd{bd) 
For x > 1, we have 

1 (-l) k 
cot -1 ^) = tan _1 (-) - x 



, -l 



3 + 



7T 



i ) 2 



x' f^ ) (2k + l)(x 2k + 1 )' 



By truncating the Taylor series after the first and second term, we have - > 
cot _1 (a;) > i — gij. Using this inequality on each summand in (|2.2|) . substituting 
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x = —j~r for each j, and then using Faulhaber's formulas for the resulting j m terms, 

•7 + 2 



we have that for b > d 



2 ' 

> Pd(b) > 



2b ™ y ' 2b 1086 3 ' 

Suppose now that b — — — a, where a is some large constant. In that case we 
have 

( ' M)> 7b~ 108p ~ 2(d 2 - na) ~ 108(rf 2 - na)* ' 

The limit of the right hand side of (|2.3() as d increases is ^ , and for all large enough 
d the right hand side is greater than J. As each Pd{b) is decreasing in b, we see 

,2 j2 

that bd > & = — — a. As we also have — > bd for large d, we have our result. □ 

Conjecture 2.4. The root of the polynomial Md(t) with largest norm, call it 7<j, 
has |Im(7d)| maximal among the imaginary parts of all roots of degree d SNN poly- 
nomials. 

Experimental data for the roots of a large number of SNN polynomials of degree 
less than or equal to seven form the basis for this conjecture. As a corollary of 
Theorem 14.21 below. Conjecture 12.41 is true when d = 2. 

3. General Bounds and the Vertical Strip Conjecture 

In this section we are interested in the following conjecture, due to the second 
author, M. Beck, J. De Loera, J. Pfeifle, and R. Stanley. 

Conjecture 3.1. (see ^) If L P (t) is an Ehrhart polynomial of degree d, then the 
roots oti of L p{t) satisfy —d < Re(ai) < d — 1 for all i. 

Wc will refer to this as the Vertical Strip Conjecture, for obvious reasons. The 
original motivation for this claim was experimental data produced from polytopes 
of relatively low dimensions, along with the fact that all real roots of a degree d 
SNN polynomial lie in the interval [— d, d— 1] (as shown in pQ). If this conjecture is 
true, one might hope that, like the norm bounds above, it actually holds for SNN 
polynomials. One step in this direction is the following. 

Theorem 3.2. For d>2, let C\ (respectively C\) be the open pointed cone in C 
with vertex d — 1 (respectively —d), angular width "^f-, and bisecting ray (d — 1, 00) 
(respectively (—00,—d)). For any SNN polynomial f of degree d, C\ U C\ does not 
contain a root of f. 

Proof. That the real roots of such polynomials are in the interval [— d, d — 1] was 
shown in [Q. We prove that the theorem holds for complex values in C\ with 
positive imaginary part. The proof is similar for the other cases. Let t G C\. Then 
the difference between the arguments of (' + ^ _ "') and ( t+<i ^"'~ 1 ) is less than ^ for 
all j, as this is equal to the difference between the arguments of [t + d — j) and 
(t — j), both of which have argument less than ^. Therefore, the points (* + ^ - "'), 
< .7 < d, lie in a common half plane with zero on the boundary. A nonnegative, 
nonzero linear combination of such points cannot be zero. □ 

Note that the Vertical Strip Conjecture follows from Theorem 13 .21 for d = 2, as 
in this case the angular widths of our cones are n. We can extend this as follows. 
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Theorem 3.3. For any SNN polynomial of degree 3 or 4, the Vertical Strip Con- 
jecture holds. 

Proof. If d = 3 or 4, we show that for every complex number z lying outside the 
vertical strip {z : -d < Rc(z) < d — 1}, the numbers ( z+ ^~ j ), < j < d, all lie 
in a half-plane with zero on the boundary. This is tantamount to showing that 
the angles A(j) formed between the vectors z — j and z + d — j , which are the 
angle differences between successive , sum to less than it for < j < d — 1. 

The situations for the real part of z being less than d and greater than d — 1 are 
symmetric (interchanging j with d—l—j), so assume that the real part of z is greater 
than d — 1. As these are real polynomials, complex roots will occur in conjugate 
pairs, hence we may assume that all relevant vectors lie in the first quadrant and 
therefore each A(j) is increased by decreasing the real part of z. Hence, we need 
to show that the angle sum ■ A(j) is at most tt when the real part of z is d — 1. 

Let z = (d — 1) + ki. An application of the Law of Cosines to the triangle with 
vertices 0, z — j, and z + d — j yields the following formula, where r = (d — 1) — j 
(the x-coordinate of z — j) and s = (2d — 1 — j) (the x-coordinate of z + d — j): 

d 2 k 2 

COSM(J) = 1 - (fc 4 + (r 2 + s 2 )fc 2 +7 ,2 s2) 

Since d is a constant, applying the AM/GM inequality to the fraction yields 
that this quantity is minimized when k 2 — rs, and increases (thus A(j) decreases) 
monotonically in both directions as k gets further from this point. 

Now, we move on to the specific applications for d = 3 and d = 4. For d = 3, 
simply applying the above formula yields (by numerical computation) : 

{A{0),A(l),A(2)) < (0.45,0.65,1.58); 

the sum of these is less than 7r, so the relevant ( Z+ ^ _ ' J ) all lie in a half-plane. 

For d = 4, we need to consider cases. Using the monotonicity results obtained 
above, we divide into the cases k 2 < 2 and k 2 > 2. For k 2 < 2 we have: 

{A(0),A(1),A{2),A(3)) < (0.24,0.38,0.70,1.58), 

while for k 2 > 2 we have: 

(A(0),A(1),A(2),A(3)) < (0.41,0.53,0.73,1.23); 

in both cases the sum of these angles is less than 7r, and so in both cases the relevant 
( Z+ d _ "0 a S am a ^ ne m a half-plane. This completes the proof that all zeroes lie in 
the vertical strip for d = 3, 4. □ 

As an aside, the above results can be extended further to show that for d = 3 or 
4, when a potential zero is not purely real, its real part must in fact satisfy tighter 
bounds than simply being between —d and d — 1. Indeed, the locus of all possible 
zeroes off the real axis is bounded by curves given by the equations = it. 

We have shown that these curves lie entirely inside the vertical strip for d = 3,4. 
These curves will be discussed more in the final section. 

While the above results are promising, in general the vertical strip conjecture is 
not satisfied by SNN polynomials. To produce examples illustrating this, we pick 
a desired root z for which the numbers ( Z+ ^ _J ) do not lie in a complex half-plane 
for a desired degree d. We can then produce a positive linear combination of these 
numbers which is equal to zero, and the polynomial encoded by these coefficients 
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will have a root equal to z. For d large enough, we can find such z which lie outside 
the vertical strip, as the following examples demonstrate. 

Example 3.4. The polynomial 

satisfies f[z) = for z w 4.00019 + 3.00963i. 

We obtained this polynomial by noting that for z = 4 + 3i, the sum of the A[j) 
exceeds ir. Therefore, there is a SNN polynomial with root 4 + 3i, and by changing 
the coefficients slightly, we obtain a polynomial with real part strictly greater than 
4. 

However, this polynomial is not an Ehrhart polynomial. In particular, it does 
not satisfy the list of inequalities satisfied by the coefficients of Ehrhart polynomials 
given in p^. Furthermore: 

Theorem 3.5. No SNN polynomial of degree 5 with a root outside the vertical strip 
can be an Ehrhart polynomial. Therefore, the Vertical Strip Conjecture holds for 
Ehrhart polynomials (though not SNN polynomials) for d = 5. 

Proof. We give a sketch of the proof. Using methods similar to those in the proof 
of Theorem 13.31 one can show that the angle sum can only barely exceed it (it is 
bounded by 3.17.) In particular, for any d outside the vertical strip but inside the 
SNN locus, the vectors Vj := ) all lie in the same half-plane for j £ {0, . . . , 4} 

and also for j £ {1, . . . , 5}. 

This means that any positive dependence among these vectors is a sum of positive 
dependences among {vo, Vj, V5} for j £ [4]. However, for all relevant z (outside the 
vertical strip but inside the SNN locus), the magnitude of v$ turns out to be much 
smaller than the magnitude of Do and V\; in particular, in each of these dependences, 
as in Example 13.41 the coefficient of v§ must be much larger than the coefficients 
of vq and v\ (if j = 1.) 

This implies that h£ > hi + h* , which violates the inequality h§ < hi + h\ given 
in p. □ 

However, by bumping up the degree, we can find candidate Ehrhart polynomials 
that may or may not be actual Ehrhart polynomials of polytopes. 

Example 3.6. Consider the polynomial g(t) of degree 26 with 

[hi,..., h* 2& ) = (1,2, 3, 4, 6, 10, 16, 27, 43, 69, 112, 181, 293, 473, 762, 0, . . . , 0) . 

Numerical approximation produces 

26.47331467- 28.51231239i 

as a root of g[t). 

For large d, there are a large number of SNN polynomials with roots outside the 
vertical strip; any point for which the angle sum is larger than n is the root 

of some SNN polynomial. If we pick a point close to the vertical strip, we will have 
a large number of SNN polynomials which have it as a root (an SNN polynomial 
is produced by any positive linear dependence among d vectors which for large d 
have roughly evenly spaced arguments). It seems as though one of these must be 
an Ehrhart polynomial, though of course discerning the Ehrhart-ness of degree 26 
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polynomials (which come from dimension 26 polytopes) with relatively large h*- 
vectors is certainly a nontrivial task. Indeed, the following question is concrete and 
unresolved. 

Question 3.7. Is the polynomial g(t) from Examvle \H.b\ an Ehrhart polynomial? 

The answer would be interesting either way; if it is an Ehrhart polynomial, it 
is a counterexample to the Vertical Strip Conjecture, and if not, new methods will 
need to be developed to verify this. 

4. Experimental Results and Bounds in Low Dimension 

For Ehrhart polynomials of degree 2, there are very tight known restrictions on 
the location of the roots. 

Theorem 4.1. (Beck, et al, see £Q.) The roots of the Ehrhart polynomial of any 
lattice 2-polytope are contained in 



For degree 2 SNN polynomials, we can produce very similar root restrictions. 
As mentioned after the proof of Theorem 13.31 the zeros of SNN polynomials in 
fixed degree are contained in regions bounded by curves given by ^ A(j) = n. By 
analyzing this curve in detail, we obtain the following. 

Theorem 4.2. The roots of any degree 2 SNN polynomial are contained in 



Proof. Let z = x + iy. Suppose that z is not pure real; consider the diagram in 
Figure ^ which shows the relevant points z — 1, z, z + 1, z + 2. Then we need to 
determine the locus of all z such that the angles 9 and <f> are complementary. The 
situation is obviously symmetric about the real axis and about x = —1/2, so we 
assume that x > —1/2,?/ > 0. In this picture, we must have 9 > <f>: the derivative 
of tan -1 x is a decreasing function of |a;|, the uncommon parts of 9 and (f> are the 
integral of this derivative over an interval of size 1, and the one which is part of 
4> is strictly further from the y-axis (since the real part of z is greater than -1/2.) 
It then follows that 9 + <f> > ir if and only if sin 9 < sin 0, or, since both sines are 
positive, if and only if sin 2 9 < sin 2 <f>. 





Z 



- 1 



Z 



z + 1 



z + 2 




Figure 1. 
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We use the Law of Sines, which states that the area of a triangle is one-half times 
the product of two adjacent sides times the sine of the included angle, on the two 
triangles with included angles 9 and <p. These triangles both have area equal to y, 
as their height is y and their base is 2. For the triangle with angle 9, we obtain: 



; \/ {x — l) 2 + y 2 \fjx + l) 2 + y 2 sin 9, or 



Ay 2 



((x~l) 2 + y 2 ))((x + l) 2 +y 2 )- 



Similarly, we obtain 



sm 



(x 2 + y 2 )((x + 2) 2 + y 2 )' 
Comparing these formulas, we algebraically manipulate: 



Ay 2 

((x-l) 2 +y 2 ))((x + l) 2 + y 2 ) < 

((x-l) 2 +y 2 )((x + l) 2 +y 2 )) > 

x 4 - 2x 2 + 1 + y 2 (2x 2 + 2) + y 4 > 

y 2 {4x + 2) + 4x 3 + 6x 2 - 1 < 

y 2 < 



V 

(x 2 + y 2 )((x + 2) 2 + y 2 ) 

(x 2 +y 2 )((x + 2) 2 +y 2 ) 

x 4 + A(x 3 + x 2 ) + y 2 (2.T 2 + Ax + 4) + y 4 



1 - Ax 3 - 6x 2 
Ax + 2 



where the last step is because Ax + 2 is positive. This reduces to y 2 < —x 2 — x+l/2, 
which is the equation of the indicated circle. 

Finally, if z is pure real, the result follows immediately from Theorem 13.21 □ 

There are also very tight restrictions on the location of roots of Ehrhart polyno- 
mials of degree 3, as the following theorem shows. 

Theorem 4.3. (Bey, et al, see [3 .) The roots of the Ehrhart polynomial of any 
lattice 3-polytope are contained in 

[-3,1] U {x + iy £ C : -1 < x < l,x 2 + y 2 < 3} . 

The curve = 7r for degree 3 SNN polynomials is shown in Figure [21 In 

degrees 2 and 3, note that the known restrictions for Ehrhart polynomial roots and 
SNN polynomial roots differ primarily in the restrictions on their real parts. 

Unfortunately, for higher degrees the bounding curves given by A(j) = tt are 
more complicated. Figure [21 shows this curve for degree 7 SNN polynomials. It is 
interesting to compare this to Figure[Sl a plot of the roots of 1000 random degree 7 
SNN polynomials, and to Figure 0] an approximation of the region containing the 
roots of degree 7 SNN polynomials as determined by Theorems 12.11 and 13.21 All 
three of these pictures have roughly the same "shape," though the random root plot 
is contained in a much smaller region than that bounded by the curve A(j) = n, 
which is a smaller region than that given by Theorems 12.11 and 13.21 
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Figure 2. The curve given by ^2,A(j) = it for degree 3 SNN polynomials. 




Figure 3. The curve given by ^A{j) = it for degree 7 SNN polynomials. 



As we have seen, both Ehrhart and SNN polynomials have roots growing quadrat- 
ically in norm as the degree grows, with the maximal roots of the extremal can- 
didates for each class differing in norm by roughly a factor of two. It would be 
interesting to see if there continues to be significant differences in the restrictions 
on the real parts of these roots in higher dimensions and if those differences can be 
quantitatively analyzed in a similar fashion. 
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Figure 4. The region containing the roots of degree 7 SNN poly- 
nomials from Theorems 12 . 1 1 and 13 . 21 




Figure 5. The roots of 1000 random degree 7 SNN polynomials. 
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